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We report on an improved scheme to generate Bose-Einstein condensates (BECs) and degenerate 
Fermi gases of strontium. This scheme allows us to create quantum gases with higher atom number, a 
shorter time of the experimental cycle, or deeper quantum degeneracy than before. We create a BEC 
of ^^Sr exceeding 10^ atoms, which is a 30-fold improvement over previously reported experiments. 
We increase the atom number of ^^Sr BECs to 2.5 x 10^ (a fivefold improvement), and refine the 
generation of attractively interacting ^^Sr BECs. We present a scheme to generate ^^Sr BECs with a 
cycle time of 2 s, which, to the best of our knowledge, is the shortest cycle time of BEC experiments 
ever reported. We create deeply-degenerate ^^Sr Fermi gases with T/Tp as low as 0.10(1), where the 
number of populated nuclear spin states can be set to any value between one and ten. Furthermore, 
we report on a total of five different double-degenerate Bose-Bose and Bose- Fermi mixtures. These 
studies prepare an excellent starting point for applications of strontium quantum gases anticipated 
in the near future. 

PACS numbers: 67.85.-d, 67.85.Fg, 67.85.Lm, 67.85.Pq 



I. INTRODUCTION 

Quantum degenerate gases of alkaline-earth-metal 
atoms, such as strontium, or of atoms with similar elec- 
tronic structure, such as ytterbium, are attracting con- 
siderable interest, fueled by the rich internal structure 
of these atoms. The two valence electrons give rise to 
singlet and triplet states, including metastable states, 
and narrow inter combination lines. In contrast to the 
bosonic isotopes, the fermionic isotopes possess a nuclear 
spin, which is crucial for many applications. Proposals 
demanding some or all of these properties describe the 
creation of artificial gauge fields [JJiSj, the implementa- 
tion of sub- wavelength optical lattices , the processing 
of quantum information [7], or the study of many-body 
systems with dipolar or quadrupolar interaction |9] . 
The nuclear spin of the fermionic isotopes is at the heart 
of proposals to study SU(A^) magnetism p^H2T] , to cre- 
ate non- Abelian artificial gauge fields [H [2] , to simulate 
lattice gauge theories [22], or to robustly store quantum 
information and perform quantum information process- 
ing [23-26 . Quantum gas mixtures of alkaline-earth- 
metal atoms with alkali atoms can be used as a basis 
for the production of ground-state open-shell molecules, 
such as RbSr [27l [28], which constitute a platform to- 
wards the simulation of lattice- spin models [29 , 30 . Bi- 
alkaline-earth- metal molecules, such as Sr2 |31|i32j, are 
sensitive and model-independent probes for variations of 
the electron-to-proton mass ratio ^1 [34|. Aside from 
degenerate gases, two- valence-electron atoms have been 
used for optical clocks [35] and other precision experi- 
ments, as well as for the production of ultracold plasmas 
[36^ . and Rydberg gases [37 . 

The first alkaline-earth-metal-like element to be cooled 
to quantum degeneracy was ytterbium [38H43j , followed 
by calcium [44 and strontium [45] [46]. Strontium 



has one fermionic and three bosonic stable isotopes, 
with vastly different natural abundances and scatter- 
ing lengths, see Tab. [ij All isotopes have been cooled to 
quantum degeneracy using slightly different approaches 
to accommodate their respective properties. The first 
isotope brought to quantum degeneracy was the low 
abundance isotope ^^Sr, since its interaction properties 
are ideal for evaporative cooling [45l|46]. Next, the ^^Sr 
isotope, which has a very small and negative scattering 
length, was Bose condensed by sympathetic cooling with 
^^Sr [47]. The last remaining bosonic isotope, ^^Sr, has a 
very large scattering length, which required evaporative 
cooling at relatively low densities to avoid three-body 
decay [48 . The fermionic isotope ^^Sr was cooled to 
quantum degeneracy in both a mixture of spin states 
[49] and as a spin-polarized sample together with ^^Sr 

In this article, we describe our scheme to create quan- 
tum degenerate samples of strontium, which improves 
several key properties substantially compared to previ- 
ously reported results. In particular, we create quantum 
gases with much larger atom number, or much shorter 
time of the experimental cycle, which increases signal 



Table I: Important properties of the four stable strontium iso- 
topes. The scattering lengths a are averages of values taken 
from [51 and 52 . Only the fermionic ^^Sr isotope has a 
nuclear spin /. 
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to noise and data rate in many applications and gives 
good prospects for using Sr as a coolant for other ele- 
ments. The temperature of degenerate Fermi gases has 
been lowered to T/Tp = 0.10(1), which is crucial for 
the exploration of SU(A^) magnetism. Furthermore, we 
report on the production of novel isotopic quantum gas 
mixtures. These improvements are obtained with a uni- 
versal experimental setup that can be used for all iso- 
topes, as well as mixtures among them, with only little 
modifications. 

The article is structured as follows: in Sec. [llj we pre- 
pare the stage by describing the experimental setup. In 
Sec. Ill we present the creation of a large ^^Sr BEC con- 
taining more than 10^ atoms and in Sec. IV a route to 
BECs containing 10^ atoms with a 2 s cycle time. In 
Sec. [V| we show improved experiments on BECs of the 
^^Sr isotope and in Sec. VI a robust route to attain BECs 
of the ^^Sr isotope. In Sec. VII we introduce two differ- 



ent double-degenerate Bose-Bose mixtures. In Sec. |VIII 
we turn to the fermionic isotope and present deeply de- 
generate samples with different numbers of occupied spin 



states. In Sec. IX we present mixtures of the fermionic 
with each of the bosonic isotopes. A summary of all 



experiments is given in Tab. Ill and a short conclusion 
follows in Sec. [Xj 



II. EXPERIMENTAL SETUP AND 
PROCEDURES 
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Figure 1: (Color online) Schematic illustration of the en- 
ergy levels and transitions used for cooling and trapping of 
Sr atoms. The blue MOT is operated on the strong ^Sq-'^Pi 
transition. The loading of the magnetic trap proceeds via 
the weak leak of the excited state (branching ratio 1:50 000) 
into the state, which itself decays with a 1:3 probability 
into the met ast able state. Here the atoms can be mag- 
netically trapped and accumulated for a long time, forming 
a reservoir. The transition ^P2-^^2 allows us to depopulate 
the metastable state by transferring the atoms into the ^Pi 
state. The latter represents the excited state of the ^Sq-^Pi 
intercombination line, used for narrow-line cooling in the red 
MOT. 



In this section, we will describe crucial parts of the 
experimental setup and give an account of the experi- 
mental sequence. The basic scheme to produce Sr quan- 
tum gases is the preparation of a laser cooled sample in a 
magneto-optical trap (MOT), the transfer of this sample 
into an optical dipole trap, followed by evaporative cool- 
ing. Laser cooling of Sr atoms has already been described 
in detail elsewhere |50l ISS'-SG and will be reviewed 
only briefly in Sec. |II A[ In Sec. |IIB| we describe the de- 
sign of a flexible optical dipole trap, which has proven 
to be the crucial element for efficient evaporation. The 
transfer of atoms from the MOT into the dipole trap is 
described in Sec. |II C| Evaporative cooling to quantum 



degeneracy is discussed in Sec. |IID| and data acquisition 
and analysis in Sec. |IIE[ A detailed description of the 
experimental setup and the laser systems can be found 
in Kefs. [571 Eg. 



A. Laser cooling 

For the preparation of a laser cooled sample of Sr in an 
ultrahigh vacuum chamber, we use a variation of stan- 
dard laser cooling and trapping techniques [59]. The 
preparation method has two stages and is based on the 
term scheme of Sr [see Fig. [l]. During the first stage, 
a transversally cooled and Zeeman-slowed atomic beam 
is captured in a "blue" MOT. The MOT cycle is not 
completely closed and pumps atoms into a metastable 



magnetic state, in which they can be accumulated in a 
magnetic trap. During the second stage, the atoms are 
pumped out of the metastable state and further cooled 
in a "red" MOT. This MOT uses a narrow-linewidth 
transition and laser cools the sample to a phase-space 
density approaching 0.1. 

All laser cooling steps of the first stage use the broad 
^So — ^Pi transition at 461 nm. The linewidth of this 
blue transition is Fbiue = 27r x 30.5 MHz, which corre- 
sponds to a comparatively high Doppler temperature of 
720 /iK. The transition is nearly closed and no repumping 
is needed to operate the blue MOT. The isotope shift of 
about 70 MHz per mass unit is larger than the transition 
linewidth, which results in isotope selective laser cool- 
ing. Operation can be switched between isotopes within 
100 ms by changing the laser frequency. Our laser source 
is a frequency doubled diode laser system with 350 mW 
output power. 

During the first stage, Sr atoms are emitted through 
a bundle of microtubes (~ 70 steel tubes, 200 /im inner 
diameter, 300 /im outer diameter, 12 mm length) from 
an oven heated to about 600° C. The resulting atomic 
beam is transversally cooled by two retroreflected laser 
beams, orthogonal to the atomic beam and to each other, 
which increases the flux of atoms captured in the MOT 
by about a factor four. The transversal cooling beams 
have an elliptical beam profile elongated along the di- 
rection of the atomic beam (15mm by 3mm waist), a 
detuning of —15 MHz from resonance, and a power of 
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Figure 2: (Color online) Schematic view of the setup as 
seen from above (a) and from the side (b). Shown are the 
quadrupole field coils (light gray areas with dark gray out- 
line), the MOT beams (rectangular blue areas), the ellip- 
tical horizontal and the circular vertical dipole trap beams 
(bow tie shaped red areas), and the direction of imaging, the 
(x — y)-direction (open arrow). The horizontal dipole trap is 
elongated in the x-direction and quantum degenerate gases 
are created in the pancake shaped crossing region of the two 
dipole trap beams. For clarity, only the focus of the hori- 
zontal (vertical) dipole trap beam is shown in the side (top) 
view. 



10 mW per beam. The atomic beam passes two consecu- 
tive differential pumping stages, helping to achieve excel- 
lent vacuum conditions in the ultrahigh vacuum cham- 
ber, which result in a lifetime of two minutes for Sr sam- 
ples stored in an optical dipole trap. The atoms are de- 
celerated by a spin-flip Zeeman-slower of 80 cm length. 
The Zeeman-slower beam has a power of ~ 35 mW, a 
detuning of —430 MHz, a waist of about 8 mm at the 
position of the MOT, and is slightly focussed onto the 
aperture of the oven. The atoms are captured by the 
blue MOT in a glass cell vacuum chamber. The MOT 
consists of three orthogonal retroreflected laser beams, 
with a detuning of -30 MHz [see Fig. [2]. The horizontal 
(vertical) MOT beams have a waist of 5 mm (4 mm) and 
a power of 4mW (ImW) each. The quadrupole mag- 
netic field of the MOT has a gradient of 55 G/cm in the 
axial direction (z-direction). 

The blue MOT cycle is not completely closed, as atoms 
in the ^Pi state can decay into the state and fur- 
ther into the metastable state, which is magnetic 
and has a natural lifetime of many minutes ^0 . Weak- 
magnetic-field seeking atoms in this state are trapped 
in the magnetic quadrupole field that is also used for 
the MOT, forming a reservoir of atoms. We operate the 
MOT until an estimated amount of 10^ atoms are cap- 
tured in this metastable state reservoir, which usually 
takes between 100 ms and a few seconds, depending on 
the natural abundance of the isotope of interest. At this 
point, all laser beams and the Zeeman slower magnetic 
field coils are switched off and the atomic beam is blocked 
by a mechanical shutter. The use of the metastable state 
reservoir for accumulation of atoms is an elegant way to 
overcome the low natural abundance of the ^^Sr isotope 
(0.56%), or to store a certain isotope while loading a 
different one. The lifetime of atoms in the reservoir is 



about 30 s. 

During the second laser cooling stage, the atoms are 
further cooled in a red MOT operated on the ^5*0 — ^Pi 
intercombination line at a wavelength of 689 nm. The 
linewidth of this transition is Fred = 27r x 7.4 kHz, 
which corresponds to a Doppler temperature of only 
180 nK. This temperature is comparable to the recoil 
temperature, which is 460 nK. The saturation intensity 
of this transition is /gat = 3/iW/cm^. The red MOT 
beams are superimposed with the blue MOT beams, 
have a waist of 3 mm, and an initial peak intensity of 
6mW/cm^ = 2000 /gat per beam. Our laser source is 
based on a extended cavity diode laser with a linewidth 
of about 2 kHz and an absolute stability better than 
500 Hz. The narrow linewidth is achieved by locking 
the laser to a cavity, which is stabilized in length using 
a spectroscopy lock on a timescale of 5 s. The required 
frequency components are derived from this master laser 
by acousto-optical modulators (AOMs) and amplified by 
slave diode lasers. 

Loading of the red MOT is accomplished by optical 
pumping on the — transition at 497 nm [see 
Fig.|l]. The repump beam has a power of a few mW and 
a waist of 8 mm. Its laser source is a frequency doubled 
diode laser, which is referenced to a stable cavity. After 
switching the repump beam on, the magnetic field gradi- 
ent is lowered to 1.15 G/cm within 0.1 ms. The tempera- 
ture of the atoms that are pumped out of the metastable 
state reservoir is on the order of the blue MOT Doppler 
temperature, which is 720 /iK. To obtain a capture ve- 
locity that is high enough for these atoms and to in- 
crease the capture volume, we frequency modulate the 
red MOT light. The modulation creates sidebands that 
cover a detuning range between —150 kHz and —8 MHz 
with a spacing of 20 kHz. The red MOT is kept in these 
conditions for 50 ms to capture the atoms. 

To increase the phase-space density, the red MOT is 
compressed by ramping the detuning, modulation width, 
and power of the MOT beams to lower absolute values 
over the next ~ 500 ms. The modulation is switched 
off already after 100 ms of this ramp. At the end of 
the ramp, the detuning of the MOT beams is typi- 
cally — 150 kHz and the peak intensity is 1.5/iW/cm^ = 
0.5 /sat per beam. This compression phase leads to a 
colder and denser sample. A few 10^ atoms can be cooled 
to temperatures of about 400 nK for the almost non- 
interacting ^^Sr isotope or about 800 nK for the other 
isotopes, including the fermionic one. The attainable 
temperature highly depends on the density of the sam- 
ple, which creates a trade-off between temperature and 
atom number. 

The operation of the fermionic red MOT is slightly 
complicated in comparison to the bosonic red MOT, be- 
cause of the existence of hyperfine structure in fermionic 
^^Sr, which is absent in bosonic Sr isotopes. Since ^^Sr 
possesses a nuclear spin of / = 9/2, its ^Pi state is split 
into three hyperfine states with F' = 11/2, 9/2, and 7/2. 
A long-lived MOT is obtained by using the transitions 
to the ^Pi{F' = 11/2) and the ^Pi{F' = 9/2) states 



4 



■ 




88Sr 








|g 










200 pm 



(b) 






200 pm 



Figure 3: (Color online) Narrowline MOTs of (a) the bosonic 
^^Sr and (b) the fermionic ^^Sr isotopes, shown by in-situ 
absorption images taken along the horizontal direction. In 
case (a), the MOT beams have a detuning of -50kHz and a 
peak intensity of /sat- Gravity and laser cooling forces bal- 
ance each other on the surface of an ellipsoid, which has a 
vertical radius of 200 /im, giving rise to a pancake shaped 
MOT. In the fermionic case (b), we operate at a detuning of 
about -20 kHz and an intensity equal to hat for both MOT 
frequency components. The atoms occupy the volume of 
an ellipsoid. In both cases, the magnetic field gradient is 
dB/dz = 1.15 G/cm, the atom number 1.3(1) x 10^, and the 
temperature ^ 700 nK. The white ellipses are a guide to the 
eye. 



simultaneously [54| l56l. During the capture and com- 
pression phases of the red MOT, both laser components 
are ramped as described above. 

After the compression phase, the shapes of the bosonic 
and fermionic red MOTs are very different [see Fig. [s]. 
In the bosonic case, the atoms occupy the surface of 
an ellipsoid, defined by the resonance condition —hA = 
^j/i^5(r), where h is the Planck constant, A the de- 
tuning of the MOT laser from the zero magnetic field 
resonance, gj = 1.5 the Lande ^-factor of the ^Pi state, 
/j^B the Bohr magneton, and B{t) the magnitude of the 
MOT quadrupole magnetic field at location r. On the 
lower part of this surface the force of gravity is com- 
parable to the restoring light force, and atoms accumu- 
late there [53, 61 . In our typical operation conditions 
(A = —150 kHz, dB/dz = 1.15 G/cm) the atoms form a 
disc with about 500 jam diameter, located 600 jam below 
the quadrupole center. By contrast, the fermionic MOT 
occupies the entire volume of an ellipsoid. In order to 
achieve densities similar to the bosonic case, we reduce 
the detuning of both MOT frequency components to a 
few linewitdhs, which reduces the size of the ellipsoid. 
The positions of the bosonic and fermionic red MOTs 
are thus different, which needs to be taken into account 
for dipole trap loading. 

Our laser cooling scheme can be easily extended to 
prepare a mixture of two Sr isotopes. For this task, the 
first cooling stage is executed twice, filling the long-lived 
metastable state reservoir consecutively with the first, 
then with the second isotope. The only difference be- 
tween these two executions is the wavelength of the blue 
laser system, which is changed from one isotope to the 
other. The second laser cooling stage is executed simul- 
taneously for both isotopes. Here, the red MOT beams 



and the repump laser contain several frequency compo- 
nents each, tuned to the transitions required for repump- 
ing and laser cooling of the two isotopes. The required 
frequency components can be easily derived from master 
lasers using AOMs. 



B. Dipole trap setup 

For further storage and cooling of the atomic cloud, we 
use a crossed beam optical dipole trap [62^ , consisting of 
a horizontal and a (near-) vertical infrared laser beam 
[see Fig.|2]. The dipole trap has its cross at the location 
of the atomic cloud and is switched on at the beginning 
of the red MOT phase. Each laser beam is derived from 
a separate multimode 5-W Yb-doped fiber laser (IPG 
YLD-5-1064-LP) operating at 1065 nm. The intensity 
of each laser beam is set, but not actively stabilized, 
by AOMs. Polarization-maintaining single-mode optical 
fibers behind the AOMs are used for mode-cleaning. 

The transfer of atoms from the red MOT into the 
dipole trap is optimized by spatial mode-matching. Since 
the bosonic red MOT is pancake-shaped, we give the hor- 
izontal dipole trap an elliptical focus, elongated in the 
horizontal plane. To shape the focus, we use a cylin- 
drical telescope between the output collimator of the 
mode-cleaning fiber and the lens focussing the beam 
onto the atoms. The beam has a vertical waist of 
18(2) /im and a horizontal waist of 250(50) /im, thus an 
aspect ratio of 15. With a typical power of 2W used 
for loading of the dipole trap, the trap frequencies are 
= 40(4) Hz, fy = 6(1) Hz, and = 620(50) Hz. The 
trap is /cb X 10(1) /iK deep, where /c^ is the Boltzmann 
constant and we take into account gravitational sagging. 

The vertical dipole trap beam is used to obtain a 
stronger confinement in the horizontal plane, but it does 
not support the atoms against gravity. This beam is fo- 
cused into the vacuum chamber by a spherical lens and 
is centered onto the focus of the horizontal dipole trap 
beam. The 1/e^-radius of the vertical beam in the plane 
of the horizontal dipole trap can be set to any value be- 
tween 20 jam and 300 /im, and is optimized separately for 
evaporative cooling of each Sr isotope or isotopic mix- 
ture. Two methods are used to adjust the beam size. 
Coarse tuning is done by changing the size of the colli- 
mated beam after the mode-cleaning fiber by exchanging 
the fiber collimator. Fine tuning is done by moving the 
position of the focussing lens along the beam axis us- 
ing a translation stage. Due to optical restrictions, this 
dipole trap beam is not exactly vertical but propagates 
at an angle of 22° with respect to gravity. To be able to 
quickly switch between two different waists, two vertical 
dipole trap beams are implemented, reaching the atoms 
under angles that differ by less than 30 mrad. In all ex- 
periments presented here, only either one or the other of 
the two beams is used. 

The role of the horizontal and vertical dipole trap 
beams for evaporative cooling is clearly distinct. The 
horizontal dipole trap beam sets the trap depth, since 
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atoms leave the trap predominantly vertically down- 
wards, aided by gravity. It also determines the verti- 
cal trap oscillation period and thereby the timescale on 
which atoms can escape the trap. Evaporative cooling 
can be limited by this timescale and therefore a high 
vertical trap frequency is desirable. The vertical dipole 
trap adds confinement in the horizontal plane and is a 
means to tune the density of the sample. This additional 
confinement is used in all experiments presented, but for 
the creation of a ^^Sr BEC with high atom number. To 
achieve a sufficiently high density for evaporative cool- 
ing, we never need horizontal trap frequencies exceeding 
the vertical trap frequency. The trap has always a pan- 
cake shape, which is fortuitous for two reasons. First, 
this trap shape is matched to the bosonic red MOT, 
which is a requirement for efficient trap loading. And 
second, evaporation happens across a large surface of 
the trap, which is good for evaporation efficiency. 

To determine the potential of the dipole trap, we per- 
form trap oscillation measurements. The measurements 
are typically performed using a nearly pure BEC in a 
trap that is slightly recompressed from the conditions in 
which the BEC was created. To measure the trap oscilla- 
tion frequencies in the radial plane, we use the additional 
vertical dipole trap beam to excite sloshing modes in the 
horizontal plane. To measure the vertical trap oscillation 
frequency, we exploit gravitational sagging and excite the 
vertical sloshing mode by a sudden compression of the 
horizontal dipole trap. From the trap frequencies and the 
power of the dipole trap lasers, we derive the 1/e^-radii 
of the dipole trap beams at the location of the cross and 
the depth of the potential, taking into account the effect 
of gravitational sagging. Knowing the radii of the dipole 
trap beams, we can calculate the trap potential for any 
isotope and dipole beam power. For cold atomic clouds 
near quantum degeneracy, the potential is usually well 
approximated by a harmonic potential. This approxi- 
mation is not always valid for atomic clouds with MOT 
temperatures, since these clouds can be larger in the hor- 
izontal direction than the size of the vertical dipole trap 
beam. 



C. Dipole trap loading and performance 

In the following, we will discuss loading of atoms into 
the dipole trap and then address two important details. 
The first detail are the frequency shifts induced by the 
dipole trap and by the magnetic field inhomogeneity and 
the second detail are the requirements for reliable dipole 
trap loading. 

To load atoms into the dipole trap, the spatial overlap 
between the red MOT and the dipole trap is fine-tuned 
by shifting the MOT quadrupole magnetic field center 
using an offset magnetic field. About 50% of the atoms 
are transferred from the MOT into the dipole trap. After 
the atoms are in the dipole trap, the MOT light is kept 
on for another 100 ms at an intensity of ~ 0.5 /gat and a 
detuning of —20 kHz from the Zeeman- and light-shifted 



cooling transition in the trap center. During this time, 
the atoms that are initially spread out in the horizontal 
dipole trap beam are pushed into the center of the trap 
by the horizontal MOT beams, thereby increasing the 
density of the sample. Afterwards all near-resonant laser 
fields are ramped off. Since the location of the bosonic 
and fermionic red MOTs are different [see Fig. [s] , prepa- 
ration of a Bose- Fermi mixture in the dipole trap requires 
a sequential loading scheme. We first load the fermions 
and then shift the quadrupole center upwards by apply- 
ing an offset magnetic field to load the bosons. We can 
transfer up to 4 x 10^ atoms into the dipole trap with 
a shot-to-shot variation of about 1%. The temperature 
of the cloud is between 400 nK and 1.5 /iK, depending 
on isotope and density. Since the ground state is non- 
magnetic for bosonic Sr, the quadrupole magnetic field 
used for the MOT does not create a potential on the 
atoms and is typically left on during evaporation. If 
fermions are used, the quadrupole field is switched off 
and a homogeneous magnetic field of 3 G is applied dur- 
ing evaporation. After loading the dipole trap, the gas is 
left to thermalize for typically 250 ms. Remarkably, the 
phase-space density at this point in the experimental se- 
quence can reach up to 0.3, which shows the power of 
narrow-line laser cooling. The lifetime of a dilute sam- 
ple, measured after some evaporative cooling and a slight 
recompression of the dipole trap, is two minutes, much 
longer than typical evaporative cooling timescales. The 
heating rate is below InK/s. 

An important effect inffuencing in-trap cooling is 
the differential light shift between the two states con- 
nected by the laser cooling transition. The different 
ac-polarizabilities of the two states together with the 
position-dependent intensity of the dipole trap beams 
leads to a position-dependent laser cooling transition fre- 
quency across the sample. For optimal laser cooling the 
differential light shift needs to be minimized. The magni- 
tude of the differential light shift depends on the dipole 
trap wavelength, the orientation of the magnetic field, 
and the polarizations of dipole trap and cooling light 
[56] . We minimize the light shift for laser cooling of a 
bosonic red MOT. In this case, the magnetic field across 
the atomic cloud is nearly homogeneous and oriented 
vertically, since the quadrupole field center is located 
600 /im above the cloud. Only the cf~ transition is near 
resonant with the cooling light. In this simple situa- 
tion and for our dipole trap wavelength of 1065 nm, it is 
possible to reduce the differential light shift to a small 
value by choosing optimal polarizations for the dipole 
trap beams. The vertical dipole trap beam has to be po- 
larized to the opposite circular direction of the vertically 
upwards propagating MOT beam. The horizontal dipole 
trap beam has to be linearly polarized in the vertical di- 
rection. 

The differential light shift is determined to within 
1 kHz using absorption imaging on the laser cooling tran- 
sition [63] and found to be positive and below +20 kHz 
for all dipole trap settings used. Since this value is larger 
than the transition linewidth of 7.4 kHz, the light shift 
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needs to be taken into account when choosing the opti- 
mum frequency for laser coohng of a sample confined in 
the dipole trap. Inconveniently, the optimum laser fre- 
quency would depend on the location of an atom in the 
trap, if the light shift across the sample would be too 
inhomogeneous. Fortunately the inhomogeneity of the 
light shift is only ~ 2 kHz, since the temperature of the 
gas during the dipole trap loading process corresponds 
to only one tenth of the trap depth. This inhomogeneity 
is smaller than the transition linewidth and comparable 
to the linewidth of our cooling laser. Therefore the laser 
cooling process at different locations in the sample is not 
significantly changed by the light shift. 

An additional frequency shift for atoms in the dipole 
trap is created by the different absolute magnetic field 
values at different locations in the trap. For the con- 
ditions under which we transfer bosonic atoms from 
the MOT into the dipole trap, this shift corresponds 
to 10 kHz across the sample, which corresponds to 1.5 
linewidths and is tolerable. 



Reliable transfer of atoms from the bosonic red MOT 
into the dipole trap requires high relative position stabil- 
ity between these two traps. This requirement is espe- 
cially strict in the vertical direction, since the horizontal 
dipole trap beam has a vertical waist of only 18 /im and 
the extension of the red MOT in the vertical direction is 
not much bigger. While the dipole trap position is very 
stable, the vertical position of the red MOT depends on 
the cooling laser frequency and on the position of the 
quadrupole magnetic field center, which both can easily 
fiuctuate. To illustrate the effect, we consider the de- 
pendence of the vertical red MOT position in absence of 
the optical dipole trap on these two parameters. With 
the vertical magnetic field gradient of 1.15 G/cm, the 
red MOT is displaced 18 /im by a change in the offset 
magnetic field of only 2 mG. The same displacement is 
reached for a frequency change of 4 kHz (see the reso- 
nance condition describing the shape of the bosonic red 
MOT above). 

We now quantify the stability of the vertical offset 
magnetic field and the laser frequency reached in the ex- 
periment. We first measure the root mean square tem- 
perature variation of a dipole trapped sample over a few 
dozen experimental runs and obtain ~ 10 nK. Then we 
measure the relation of the temperature with the de- 
tuning of the cooling light and obtain 25nK/kHz over a 
range of 20 kHz. Assuming that only the laser frequency 
fiuctuates, we deduce an absolute laser frequency stabil- 
ity of better than 400 Hz. The same temperature vari- 
ation can be caused by a vertical offset magnetic field 
fiuctuation of 200 /iG. Because of the light shift induced 
by the dipole trap on the cooling transition, also a slight 
change of the polarization of the dipole trap beams could 
contribute to the observed temperature variation. 



D. Evaporative cooling 

To increase the phase-space density and achieve quan- 
tum degeneracy, we perform forced evaporative cooling. 
The intensities of the dipole trap beams are lowered in 
an approximately exponential manner. The intensity of 
the horizontal dipole trap beam sets the trap depth and 
the intensity of the vertical dipole trap beam is used to 
tune the density of the sample. Evaporative cooling to 
quantum degeneracy takes between 550 ms and 26 s. 

We now discuss some of the conditions required for 
successful evaporative cooling and the parameters with 
which we can infiuence these conditions. Evaporative 
cooling requires a high rate of elastic collisions to pro- 
ceed faster than competing heating processes and loss 
of atoms. An important contribution to atom loss are 
three-body collisions. Both, elastic collisions and three- 
body collisions, depend on the scattering length a and 
the density n. The different Sr isotopes and isotope mix- 
tures span a large range of elastic scattering lengths [see 
Tab. [H]. These scattering lengths are not tunable in a 
way suitable for evaporative cooling, since magnetic Fes- 
hbach resonances do not exist and optical Feshbach res- 
onances were shown to be accompanied by strong losses 
j64H67] . This situation leaves the density as the only 
practical parameter to increase the elastic collision rate 
and achieve a high ratio of elastic to three-body colli- 
sions. The elastic scattering rate scales as a^n, which 
favors large densities. The three-body collision rate has 
an upper limit proportional to a^v? [GS*, |69]. Therefore 
the ratio of elastic to three-body collisions is approxi- 
mately l/{o?n^) and favors low densities, especially for 
the isotope with large scattering length, ^^Sr. A limit 
to lowering the density is set by other loss processes and 
heating. Importantly, the optimum density depends also 
on the scattering properties of each isotope or isotopic 
mixture and is therefore tuned experimentally by adapt- 
ing the power and the radius of the vertical dipole trap 
beam. 



E. Data acquisition and analysis 

We characterize the gas by analysis of absorption im- 
ages [see e.g. Fig.ga)]. These images are taken in the 
(x — y)-direction [see Fig. [2] after a free expansion time 
between 20 ms and 30 ms using the ^Sq — ^Pi transition. 
For high-density samples we detune the absorption imag- 
ing beam to reduce the optical density. The quadrupole 
magnetic field used for the MOT is switched off at the be- 
ginning of the free expansion. To analyze the images, we 
apply two-dimensional fits. For bosons, during the early 
stage of evaporation, the momentum distribution of the 
atom cloud is thermal and can in most cases be described 
by a single Gaussian. After the BEG phase-transition 
has occurred, we employ a bimodal fit to our data, con- 
sisting of a Thomas-Fermi distribution describing the 
BEG and a Gaussian for the thermal atoms [70 . The 
temperature of the sample is extracted from the Gaus- 
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sian part of the fit. For fermions we use a Fermi-Dirac 
distribution, as in our previous work [50 . When work- 
ing with a mixture of two isotopes, we only detect either 
one or the other of the isotopes after each experimen- 
tal cycle. To avoid cross talk between the two isotopes 
from off-resonant absorption, we remove the isotope we 
do not want to image after 17 ms of free expansion by 
a blast of resonant light on the very isotope-selective 
-'^S'o — ^Pi inter combination line. To detect the inter- 
nal state distribution of fermionic samples, we use the 
optical Stern- Gerlach technique [43, 63 . For situations 
in which the dipole trap potential can be well approxi- 
mated by a harmonic potential, we calculate quantities 
of interest, as for example the peak density or the elas- 
tic collision rate, from trap oscillation frequencies, atom 
number, and temperature [7QH73] . 



III. 30-FOLD INCREASE IN ^^SR BEC ATOM 
NUMBER 

The first Sr isotope cooled to quantum degeneracy was 
^^Sr and BECs of 3 x 10^ atoms have been obtained 
[45] |46] . Here we increase this atom number by a fac- 
tor 30, creating BECs exceeding 10^ atoms. To overcome 
the low natural abundance of ^^Sr, we accumulate atoms 
in the metastable reservoir for 40 s. This time is slightly 
longer than the lifetime of the gas in the reservoir, and 
further loading does not increase the atom number sig- 
nificantly. The atoms are returned into the ground state, 
cooled and compressed by the red MOT, and transferred 
into the dipole trap. For this experiment, we use only the 
horizontal dipole trap beam, which has an initial depth 
of ks X 12(1) fiK and provides initial trapping frequen- 
cies of = 45(5) Hz, fy = 6(1) Hz, and f, = 650(50) Hz. 
After ramping the red MOT light off over 100 ms, the gas 
is allowed to thermalize in the dipole trap for 250 ms. At 
this point, about 4 x 10^ atoms are in the dipole trap at 
a temperature of 1.5 /iK. The peak density of the gas is 
7 X lO^^cm"^, the average elastic collision rate 650 s~^, 
and the peak phase-space density 0.3. The power of the 
dipole trap is reduced exponentially from its initial value 
of 2.4 W to 425 mW within 10 s. 

We investigate the behavior of the gas during evap- 
orative cooling by careful analysis of absorption images 
taken after a free expansion time of 30 ms [see Fig. [4]. To 
reduce the optical density the imaging light is detuned 
by 48 MHz from the - ^Pi transition. Figure [i^b) 
shows the evolution of thermal and BEC atom number 
and Fig. Qc) the evolution of temperature during evap- 
oration. After 7 s of evaporation a BEC is detected. At 
this time, 2.5 x 10^ atoms remain in the trap at a tem- 
perature of about 400 nK. The evaporation efficiency is 
high with four orders of magnitude gain in phase-space- 
density for a factor ten of atoms lost. After 10 s of evap- 
oration, we obtain an almost pure BEC of 1.1(1) x 10^ 
atoms. The trap oscillation frequencies at this time are 
= 20(3) Hz, fy = 2.5(5) Hz, and f, = 260(50) Hz. 
The BEC has a peak density of 2.2 x 10^^ cm~^ and the 




(b) 

40x 10^ 




-500 500 

horizontal position ([jm) 

(c) 

1.51- 



2 4 

evaporation time (s) 



2 4 
evaporation time (s) 



8 10 



Figure 4: (Color online) Generation of ^^Sr BECs exceed- 
ing 10^ atoms. Panel (a) shows a series of absorption images 
(left) and density profiles (right) for different times along the 
evaporative cooling ramp. The density profiles are obtained 
from the absorption images by integration along the vertical 
direction. The black circles are data points, which are fitted 
with a bimodal distribution (solid red line), capturing the 
thermal fraction (dashed blue line) and the BEC. Panel (b) 
shows the evolution of the atom number in the thermal (black 
squares) and condensate fraction (red circles) during evapo- 
ration. The temperature of the thermal fraction is shown in 
(c). 



shape of an elongated pancake with Thomas-Fermi radii 
of about Rx = 40 /im, Ry = 300 /im, and Rz = 3 /im. 
During the time of flight, the BEC expands by far the 
most in the vertical direction. Under the direction of 
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imaging used, the shape of the resulting elUpsoid ap- 
pears nearly circular [see Fig. |4ja)]. The lifetime of the 
BEC is 15 s, likely limited by three-body loss. 

The atom number of this BEC is a factor 60 larger 
than in our previously reported work [45], and is to 
the best of our knowledge the largest BEC ever created 
by evaporative cooling in an optical dipole trap. The 
changes to our apparatus that led to this improvement 
are a higher oven flux, a reduction of the linewidth and 
an increase in the frequency accuracy of the red MOT 
cooling laser, an improved dipole trap design, a new laser 
source for the dipole trap with less intensity noise, and 
a careful optimization of the experimental sequence. It 
is remarkable that the Sr BEC with the largest atom 
number is created with the Sr isotope of lowest natural 
abundance. 

An increase of the BEC atom number to nearly 10^ 
should be achievable by simple improvements. With 
more power available for the horizontal dipole trap, its 
horizontal waist could be increased while keeping the 
trap depth constant. The volume of the dipole trap 
would thus be increased and the trap oscillation fre- 
quency in the x-direction decreased. To obtain the same 
initial peak density, more atoms would have to be loaded 
into the dipole trap. This task could be achieved by 
improving the oven design to increase the atomic flux. 
An increased density of the final red MOT needs to be 
avoided to limit atom loss by light-assisted collisions, 
which could be done by using a larger detuning of the 
cooling laser, which increases the radial size of the red 
MOT [see Fig. [s]. A limit of this strategy to increase 
the atom number occurs when the trap oscillation fre- 
quency in the x-direction becomes too low at the end of 
evaporation to permit the formation of a phase-coherent 
BEC during its lifetime. A decrease to the value cur- 
rently used in the y-direction (2.5 Hz) is certainly pos- 
sible, which would already admit to increase the BEC 
atom number by a factor of six. 

The large Sr BEC created here will increase the sig- 
nal to noise in any study using ^^Sr and is a promising 
coolant for other isotopes or atomic species. 



IV. GENERATION OF BEGS WITH A GYGLE 
TIME OF TWO SECONDS 

In the previous section, we have reported on experi- 
ments optimized for a large number of atoms in the BEC. 
We can also optimize our experimental sequence for a 
short cycle time. Nearly all experiments profit from the 
higher data rate made possible by a shorter cycle time. 
Precision measurement devices, such as atom interfer- 
ometers, do require high repetition rates or a favorable 
ratio of probe time versus cycle time and might profit 
from the coherence of a BEC. Quantum gas experiments 
taking place in an environment of poor vacuum quality 
also require a short production time. Most experiments 
have cycle times of a few ten seconds. Experiments that 
have been optimized for speed while using an all-optical 
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Figure 5: (Color online) Generation of ^^Sr BECs with a 
short cycle time of the experimental sequence. Panel (a) 
shows a sketch of the experimental sequence. Absorption im- 
ages taken after 22 ms free expansion and integrated density 
profiles show the appearance of an an almost pure BEC of 
10^ atoms at the end of evaporation (b). Dashed blue lines 
indicate Gaussians, dotted green lines double- Gaussians, and 
solid red lines bi- and trimodal distributions, see the text for 
details. Panel (c) shows the development of the total atom 
number and of the temperature corresponding to the nar- 
rower of the two Gaussian distributions during evaporation. 
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approach achieve cycle times of 3 s for degenerate bosonic 
[44^,77 and lis for fermionic quantum gases |75j. The 
same cycle time can be reached by using magnetic trap- 
ping near the surface of a microchip [T^ [77] . 

Making use of the very high phase-space density 
achieved already in the red MOT, as well as the excel- 
lent scattering properties of ^^Sr, we are able to reduce 
the cycle time to 2 s [see Fig. [sja)]. At the beginning 
of the cycle, we operate the blue MOT for 800 ms to 
load the metastable reservoir. A short flash of repump 
light returns the metastable atoms into the ground state, 
where they are trapped, compressed, and cooled to about 
1.2 /iK by the red MOT. Close to 4 x 10^ atoms are 
loaded into a dipole trap, which is formed by the hor- 
izontal sheet and a vertical beam of 25 jam 1/e^-radius 
in the plane of the horizontal dipole trap. The atomic 
cloud is not only populating the cross of the dipole trap, 
but extends ~ 1 mm along the horizontal dipole trap. 
Forced evaporation reduces the trap depth over 550 ms 
with an exponential time constant of about 250 ms. 

Absorption images are taken after 22 ms free expan- 
sion and analyzed as described in Sec. HE [see Fig.jsjb)]. 
In this experiment we find that a single Gaussian is not 
sufficient to describe the thermal fraction. After about 
100 ms of evaporation, the distribution begins to deviate 
from a Gaussian without being close to quantum degen- 
eracy [see Fig. [sjb, case i)], indicating that part of the 
atoms reach a lower temperature than the rest of the 
sample. These are the atoms located in the crossing 
region of the horizontal and vertical beams, where the 
density is higher and evaporation is more efficient. After 
further evaporation, two samples with different tempera- 
tures are clearly visible, and we use a double Gaussian to 
fit our data (case ii). The phase transition occurs after 
about 270 ms of evaporation, and we employ a trimodal 
fit to our data, describing the BEG, the thermal fraction 
of atoms in the crossing region, and the thermal atoms 
in the horizontal beam (cases iii and iv). After about 
480 ms of evaporation, the thermal fraction within the 
crossing region cannot be discerned, indicating an es- 
sentially pure BEG in this region. For this last stage 
of evaporation, we employ a bimodal fit to the data, 
capturing the BEG and the thermal atoms outside the 
crossing region (case v). Further evaporation does not 
increase the BEG atom number, but efficiently removes 
the thermal atoms. The evolution of atom number and 
temperature during evaporation is shown in Fig. [sjc). 
About 10^ atoms reside in the BEG at the end of evap- 
oration. 

The read-out of the charged-coupled device (GGD) 
chip used for imaging can be performed during the con- 
secutive experimental cycle and is therefore not included 
in the 2 s period. The cycle time could be improved sub- 
stantially if the reservoir loading time (800 ms in this 
experiment) was reduced, e.g. by increasing the oven 
ffux. It seems that cycle times approaching 1 s are within 
reach. 



V. FIVEFOLD INCREASE IN ^^SR BEC ATOM 
NUMBER 

Some isotopes of alkaline-earth atoms feature large 
positive scattering lengths, such as ^^Ca, ^^Ca, "^^Ca [78], 
and ^^Sr. While scattering between atoms provides ther- 
malization during evaporation, there is a downside of a 
very large scattering length a: Inelastic three-body losses 
have an upper limit proportional to |68l [69] , and can 
reduce the evaporation efficiency drastically. Magnetic 
Feshbach resonances, a widely used means to tune the 
scattering length in ultracold samples, are absent in the 
alkaline-earth species, and a different strategy to reach 
degeneracy despite the large scattering length is needed. 

In previous work [48 , we have shown degeneracy of 
^^Sr, which has a scattering length of about 800 ao 
[51 . In that experiment, the crucial innovation was to 
perform evaporation at a comparatively low density of 
3 X 10^^ cm~^ in a dipole trap of large volume. Two- 
body collisions, vital for thermalization, scale propor- 
tional with the density n, while detrimental three-body 
collisions scale as n^. At small enough densities, evapo- 
ration can be efficient even for large scattering lengths. 
The trap was oblate with initial trap frequencies of 
= 30 Hz, fy = 3 Hz, and = 260 Hz. After 4.8 s of 
evaporation a nearly pure BEC containing 5000 atoms 
were created. 

Improving this experiment, we are now increasing the 
vertical trap frequency by a factor two, which allows us 
to evaporate faster. Using a 500 ms blue MOT stage, we 
load 9 X 10^ atoms at a temperature of about 1 /iK into 
the dipole trap. The initial density is about 10^^ cm~^ 
and the average elastic collision rate 200 s~^. We per- 
form evaporation much faster than in the previous ^^Sr 
experiment, in just 800 ms, which helps to avoid strong 
atom loss from three-body collisions. The onset of BEG 
is observed after 600 ms of evaporation at a tempera- 
ture of about 70 nK with 350 000 atoms present. Further 
evaporation results in almost pure BECs of 25 000 atoms, 
which is an improvement of a factor five in atom number 
over our previous results [see Fig. [6]^ a)]. The cycle time 
of this experiment is again short, just 2.1s. 



VI. RELIABLE GENERATION OF AN 
ATTRACTIVE ^^SR BEC 

The isotope ^^Sr is by far the most abundant one, and 
early research on cold gases of Sr focused on this isotope. 
The 5- wave scattering length a of ^^Sr turned out to be 
very close to zero, which makes it an excellent choice 
for precision measurements |79l|80]. Optical clocks, for 
example, can suffer from density-dependent mean-field 
shifts, which scale cx a. Similarly, force sensors based on 
Bloch oscillations [80 suffer from dephasing by elastic 
collisions [81. No other atomic species that has been 
laser-cooled so far combines the properties of high nat- 
ural abundance, non-magnetic ground state, and very 
small interaction strength. 
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Figure 6: (Color online) Essentially pure BECs of (a) ^^Sr 
and (b) ^^Sr. The free expansion time is 25 ms for both im- 
ages, and the geometric and color scales are identical. The 
BEC of ^^Sr contains 25 000 atoms and expands to much 
larger size compared to the almost non-interacting BEC of 
^^Sr, which contains 5000 atoms. 



It is this small (and slightly negative) scattering length 
that hampered evaporative cooling of ^^Sr towards de- 
generacy for a long time [82j |83] . Sympathetic cooling 
with the second- most abundant isotope, ^^Sr, could not 
reach quantum degeneracy because of the large 3-body 
recombination rate of ^^Sr [8l] . 

Efficient evaporation to quantum degeneracy was first 
achieved by the Rice group |47] using the fermionic iso- 
tope ^^Sr as cooling agent. An ultracold gas of ^^Sr 
is slightly attractive for T ^ 0, and indeed Mickelson 
and co-workers observed a limited number of atoms in 
the BEC, in good agreement with a simple mean- field 
model 

Here, we confirm the work of the Rice group and high- 
light a few differences of our own experimental work. We 
start with an isotopic mixture in a crossed-beam dipole 
trap, containing 4.5 x 10^ atoms of ^^Sr at 560 nK and 
1.9 X 10^ atoms of ^'^Sr at 1.0 jaK. The sample thermal- 
izes within the first second of forced evaporative cooling, 
and remains thermalized for the remainder of the 4-s 
ramp. The BEC of ^^Sr appears after 2.1 s, with 2.7 x 10^ 
atoms present at a temperature of 320 nK. Further evap- 
oration does not increase the BEC atom number beyond 
a few thousand atoms, but the thermal fraction can be 
removed to obtain very pure BECs at the end of evapo- 
ration [see Fig. [6]^b)]. 

A homogeneous BEC with negative scattering length 
is unstable for all atom numbers, but a trapped BEC can 
be stabilized through the zero-point kinetic energy of the 
trap up to a maximum atom number of Nc = ka\^o/a 
[TT. Here, aho = \/h/{muj) is the harmonic oscillator 
length, Co = 27r{fxfyfz)^^^ and the prefactor k ~ 0.55 is 
only slightly dependent on the trapping geometry [86] . 
The maximum atom number thus depends only weakly 
on the trap frequencies, and one cannot expect substan- 
tial improvements from a change in trapping geometry. 
As the BEC atom number increases and eventually sur- 
passes this limit during evaporation, part of the atoms 
are ejected out of the BEC [87 , leaving the remnant 
behind to form a stable BEC again. This cycle contin- 



ues throughout evaporation, until the remaining atom 
number is smaller than Nc or the elastic collision rate 
is too low to reduce the temperature further and add 
atoms to the BEC. Consequently, one would expect a 
random atom number N < Nc in the BEC at the end 
of evaporation. This behavior is close to the observation 
of Ref. [47 : for a large number of experimental realiza- 
tions, the BEC atom numbers were evenly distributed 
between ~ 0.4 TVc and Nc. 

In view of the application of degenerate ^^Sr in pre- 
cision measurements, we will now show how the shot- 
to-shot variation in BEC atom number can be reduced. 
We deliberately start with a low initial number of ^^Sr 
atoms to perform inefficient evaporation. After 80% (or 
3.5 s) of evaporation, only 12 000 ^^Sr atoms remain, and 
the scattering rate of a thermal ^^Sr atom with an ^^Sr 
atom has dropped from initially 250 s~^ to below 4s~^. 
At this point, we do not expect any appreciable cooling 
to take place. Instead, further lowering of the trap depth 
removes all thermal atoms from the trap, leaving behind 
a pure BEC of 6000 atoms. The atom number of this 
BEC is certainly below Nc = 7500 for our conditions of 
uj = 27T X 50(5) Hz and a scattering length of a = — 2ao 
[52 . The shot-to-shot variation in BEC atom number, 
however, is surprisingly small: We obtain a standard de- 
viation of only 9% in the atom number. This value is 
certainly larger than our typical ^^Sr BEC atom number 
fluctuation of about 1%, but it constitutes a substantial 
improvement towards precision measurements with ^^Sr 
BECs. For a higher initial ^^Sr atom number, we never 
observe a BEC exceeding 7000 atoms. When identifying 
this number with Nc^ we can deduce a scattering length 
of — 2.2(2) ao, which agrees well with the value derived 
from two-photon spectroscopy [52] . 

In conclusion, we have confirmed the ^^Sr BEC exper- 
iments by the Rice group. We observe a much smaller 
variation in BEC atom number at the end of evapora- 
tion, which might be attributed to a slightly different 
choice of initial conditions. 



VII. DOUBLE-DEGENERATE BOSE-BOSE 
MIXTURES OF SR 

Mixtures of two Bose-degenerate gases of different iso- 
topes or elements allow the study of interesting phenom- 
ena, such as the miscibility and phase separation of two 
quantum fluids [88l[9Q] or the influence of one compo- 
nent on the other during the superfluid-to-Mott insula- 
tor transition |91]. They allow the application of iso- 
tope or species specific optical dipole potentials, which 
have been used to study entropy exchange between two 
bosonic gases [92], or the dynamics of impurities in a 
one-dimensional Bose gas [93]. Within the alkalis, a 
number of quantum degenerate Bose-Bose mixtures have 
been realized: "^^K + ^'^Rb [94 , ^^Rb + ^'^Rb [95,, and 
8^Rb + i^^Cs [96]. 

The many bosonic isotopes of alkaline-earth-metal-like 
elements in principle allow the creation of many different 
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Figure 7: (Color online) Creation of double- degenerate Bose- 
Bose mixtures. The development of atom number (a) and 
temperature (b) during evaporation are shown for the mix- 
ture of ^^Sr (black squares) and ^^Sr (red circles). Despite a 
significant temperature difference at the start of evaporation, 
the sample thermalizes within 1 s. Below, we show absorp- 
tion images of the mixtures ^^Sr + ^^Sr (c) and ^^Sr + ^^Sr 
(d), taken after 26 ms of free expansion. The BEC of ^^Sr 
is nearly pure and contains 2 x 10^ atoms. The other three 
BECs are accompanied by a significant amount of thermal 
gas. The insets show the optical density integrated along the 
vertical direction, with bimodal fits to the data. Imaging of 
^^Sr is performed with a detuning of 30 MHz, the other iso- 
topes are imaged on resonance. The color scale of the ^^Sr 
image differs from the color scale used for the other three 
images, but the integrated optical densities in the insets are 
to scale. 



the most abundant Ca isotopes are quite large [44| 178], 
In Yb, two out of five bosonic isotopes have large neg- 
ative scattering lengths [97], excluding many possible 
combinations of isotopes. One remaining combination, 
^^^Yb + ^^"^Yb, has a large and negative interspecies 
scattering length. One of the two remaining combina- 
tions, ^^^Yb + ^^^Yb, has been brought to double degen- 
eracy very recently, with 9000 atoms in the BEC of each 
species [42 . The interspecies scattering length between 
these two isotopes is 2.4(3.4) ao and provides only mi- 
nuscule interaction between the two. The three bosonic 
isotopes of Sr give rise to three different two-isotope mix- 
tures [see Tab. [ll]. Of these the mixtures, ^"^Sr + ^^Sr 
suffers from a large interspecies scattering length. 

Here we report on double-degenerate Bose-Bose mix- 
tures of the combinations ^^Sr + ^^Sr and ^^Sr + ^^Sr, 
which have interspecies scattering lengths of 32 ao and 
97 ao, respectively. The experimental realization is 
straightforward: We consecutively load the two isotopes 
into the reservoir, repump them simultaneously on their 
respective — transitions, and operate two red 
MOTs simultaneously. The mixture is loaded into the 
dipole trap and subsequently evaporated to form two 
BECs. Imaging is performed on the blue ^Sq — ^Pi tran- 
sition, and we image only one isotope per experimental 
run. The frequency shift between the isotopes is only 
about 4.5 linewidths. To avoid a contribution of the 
unwanted isotope to the absorption image, we remove 
the unwanted species by an 8-ms pulse of resonant light 
on the very isotope selective ^Sq — ^Pi intercombination 
transition. To avoid a momentum distribution change of 
the imaged species by interspecies collisions, the pulse of 
light is applied after 17 ms of free expansion, when the 
density of the sample has decreased significantly. 

The experimental results are shown in Fig. [7[ We 
will discuss the ^^Sr + ^^Sr combination first: 2.3 x 10^ 
(3.3 X 10^) atoms of ^^Sr (^^Sr) are loaded into the dipole 
trap, consisting of the horizontal beam and a weak ver- 
tical beam for additional axial confinement. The ini- 
tial temperatures of the two species are quite differ- 
ent: 950 nK for ^^Sr and 720 nK for ^^Sr, which reflects 
the different intraspecies scattering behavior. The inter- 
species scattering length is around 100 ao, and the two 
species clearly thermalize to reach equilibrium after 1 s 



Bose-Bose mixtures. Unfortunately, for many of these 
mixtures, the interaction properties are unfavorable to 
create large and stable BECs. The absolute value of the 
two intra- and the interspecies scattering length must 
not be too large to avoid rapid decay and phase sep- 
aration, but large enough for thermalization. The in- 
traspecies scattering lengths should not be strongly neg- 
ative to permit the formation of detectably large BECs 
[41]. The scattering length of alkaline-earth- metal-like 
atoms can only be tuned by optical Feshbach resonances, 
which introduce losses [65 -67J. These limitations reduce 
the number of possible binary mixtures considerably. In 
particular, all combination of bosonic Ca isotopes seem 
unfavorable, since all intraspecies scattering lengths of 



Table II: Scattering lengths a between the Sr isotopes, given 
in units of ao. The values are averages of the values given 
in [51] and [52], and the uncertainty is a few ao except for 
the two very large values (a > 500 ao), where the uncertainty 
is much larger. All mixtures of isotopes, of which double- 
degenerate samples are presented in this work, are marked in 
bold. The natural abundance is given in the last column. 
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of evaporation. As the trap depth is lowered further, 
we observe the onset of BEC in ^^Sr (^^Sr) after 2.0 s 
(2.3 s). At the end of our evaporation ramp, which lasts 
2.4 s, we obtain 10 000 (3000) atoms of ^^Sr (^^Sr) in the 
condensate fraction [see Fig. [7|^a)]. Further evaporation 
does not increase the BEC atom numbers. 

In a second experiment, we investigate the ^^Sr + ^^Sr 
mixture with an interspecies scattering length of 32 ao . 
Starting out with 10 x 10^ (1.5 x 10^) atoms of ^"^Sr (^^Sr) 
in the dipole trap, we perform forced evaporation over 
2 s, and the two species remain in perfect thermal equi- 
librium throughout this time. The phase transition of 
^'^Sr is observed already after 1.3 s, with about 2.5 x 10^ 
atoms present at a temperature of 200 nK. After 1.9 s, 
the BEC is essentially pure and contains up to 2 x 10^ 
atoms. The atom number of ^^Sr is kept considerably 
lower to avoid three-body loss. The phase transition 
occurs later: after 1.7 s, with 4 x 10^ atoms at a temper- 
ature of 130 nK. Till the end of evaporation, the BEC 
fraction grows to 8000 atoms but remains accompanied 
by a large thermal fraction [see Fig. [7|^d)]. 

In conclusion, we have presented two binary Bose- 
Bose mixtures of alkaline-earth atoms with appreciable 
interaction between the two species. These mixtures en- 
joy the property that isotope-selective optical traps can 
be operated close to one of the intercombination lines. 
This might allow for an individual addressing of the iso- 
topes by a dipole trap operated close to these transitions 
[6] , reminiscent of the case of Rb in its hyperfine states 
|F = 1) and \F = 2) [98 or nuclear substates in Yb and 
Sr[43l[63]. 



VIII. DEEPLY-DEGENERATE FERMI GASES 
OF ^^SR 

A wealth of recent proposals suggest to employ the 
fermionic isotopes of alkaline-earth-metal-like atoms as 
a platform for the simulation of SU(A/') magnetism [10- 
[21], for the generation of non-Abelian artificial gauge 
fields [Hll], to simulate lattice gauge theories [22 , or for 
quantum computation schemes [23r>26 . 

Elements with a large nuclear spin are especially well 
suited for some of these proposals. They allow to en- 
code several qubits in one atom aiid could lead to 
exotic quantum phases, as chiral spin liquids, in the con- 
text of SU(A^) magnetism [T3^. Furthermore, it has been 
shown that the temperature of a lattice gas is lower for a 
mixture containing a large number of nuclear spin states 
after loading the lattice from a bulk sample [20^, '21] [99] . 
The largest nuclear spin of any alkaline-earth-metal-like 
atom is 9/2, and it occurs in the nuclei of ^''Sr and of 
two radioactive nobelium isotopes. This fact makes ^^Sr 
with its ten spin states an exceptional candidate for the 
studies mentioned. 

The study of SU(A/') magnetism in a lattice requires 
the temperature of the sample to be below the super- 
exchange scale, /U ^ where t is the tunnel matrix ele- 
ment and U the on-site interaction energy [21^. A high 




Figure 8: (Color online) Detection of spin state distribution 
using the optical Stern-Gerlach technique 43 , 63 . Samples 
of ^^Sr in a ten- state mixture or optically pumped into two 
or one spin states are shown. 



degree of degeneracy in the bulk would constitute a good 
starting point for subsequent loading of the lattice. A 
value of T/Tp = 0.26(5) in a ten-state mixture [49^ and 
of T/Tf = 0.30(5) in a single-spin state [50] of ^^Sr 
were presented already, both with about 20 000 atoms 
in each Fermi sea. Neither of these two experiments 
had the means to set the number of spin states at will. 
Here, we present degenerate Fermi gases of Sr at consid- 
erably deeper degeneracy with an arbitrary number of 
spin states between 1 and 10. 

About 5 X 10^ atoms of ^^Sr are loaded into the optical 
dipole trap, where we measure a temperature of 1.2 /iK. 
At this time, atoms are in a roughly even mixture of all 
spins. To prepare the desired spin mixture, we perform 
optical pumping on the ^5'o - ^Pi, |F = 9/2) ^ = 
9/2) transition at a small guiding field of 3G, which 
splits adjacent mi?/ -states by 260 kHz, corresponding to 
35 linewidths. As we have shown before [63], we can 
prepare any combination and relative population of the 
ten spin states. The optical pumping is optimized using 
the optical Stern-Gerlach technique [see Fig. [s] [43l (63]. 
and quantified using state-selective absorption imaging 
on the intercombination line [63^. We can reduce the 
population of undesired spin states to below 0.1%, where 
this value is limited by our detection threshold of 3000 
atoms. After the spin preparation, which does not heat 
the sample, we perform evaporative cooling. Evapora- 
tion proceeds in two stages. The first stage lasts 16 s, 
during which the power of the horizontal dipole trap 
beam is reduced by a factor of five. At this point, the 
gas enters the degenerate regime with a typical temper- 
ature of around 0.3 and evaporative cooling becomes 
less efficient because of Pauli blocking [100 . We com- 
pensate for this effect by a reduced ramp speed during 
the second evaporation stage. The power of the hori- 
zontal beam is only reduced by a factor 1.5 during 10 s. 
Slightly different final trap depths are used for different 
numbers of populated spin states. Trap frequencies at 
the end of evaporation are fx ~ 30 Hz, fy ~ 30 Hz, and 
fz ~ 200 Hz. The case of a single populated spin state 
is unique since identical fermions do not collide at low 
temperature because of symmetry requirements, imped- 



13 




4 6 
atom number per spin state 



10x10' 



Figure 9: (Color online) Deeply degenerate Fermi gases of 
^^Sr in a balanced mixture of ten nuclear spin states. Panel 
(a) shows the azimuthally averaged density distribution of 
a degenerate Fermi gas at T/Tp — 0.08(1) after 25.4 ms of 
free expansion (black circles). The measurement is well de- 
scribed by a Fermi-Dirac distribution (red line) but not by a 
Gaussian (dashed blue line). The corresponding absorption 
image is shown in the inset. Panel (b) shows the evolution of 
the atom number per spin state and of the temperature dur- 
ing evaporation, derived from two-dimensional fits of Fermi- 
Dirac distributions to time-of-flight absorption images. Panel 
(c) shows the value of T /Tp in dependence of atom number 
per spin state, derived from the fugacity (open circles) or 
calculated from the temperature, the atom number, and the 
trap oscillation frequencies (filled circles). Some data points 
have been taken multiple times (red circles with error bars) 
to determine the statistical uncertainty. 



sorption pictures as in our previous work [50^. Two 
methods are used to determine T/Tp- Either T is 
determined by the fit and Tp is calculated from the 
atom number A^at and average trap frequency as Tp = 
hCd{6N^ty^^/kB. Or T/Tp is calculated directly from 
the fugacity, which is a fit parameter. The momentum 
distribution of a gas at small T/Tp strongly deviates 
from a Gaussian shape, which we show in azimuthally 
integrated profiles [see Fig. |9ja)]. 

We will limit the presentation of our data to the cases 
of N = 10, N = 2, and N = 1. For the ten-state mix- 
ture, the fraction of Pauli-forbidden collisions is small, 
and evaporative cooling performs well to yield a stack 
of ten spatially overlapping Fermi seas. Figures |9jb) 
show the development of atom number per spin state and 
temperature along the evaporation ramp. Figure ^c) 
shows T/Tp in dependence of atom number. With about 
30 000 atoms per spin state, we obtain T/Tjr = 0.10(1) at 
Tp = 160 nK. The errors given here are statistical errors 
of multiple experimental realizations, and we estimate 
systematic errors to be of similar magnitude. 

For the two-state mixture, we pump all atoms into a 
balanced population of mp = +9/2 and mp = +7/2 
states. The total atom number, initial temperature, and 
evaporation trajectory are identical to the previous case 
of a ten-state mixture, but there is a crucial difference: 
Only half of all collisions possible for distinguishable par- 
ticles are Pauli-allowed in the binary mixture, leading to 
a decreased thermalization rate. The reduced evapora- 
tion efficiency immediately manifests itself in the degree 
of degeneracy reached: Despite the higher atom number 
per spin state, we can reach only T/Tp = 0.20(1) with 
1.0 X 10^ atom per spin state and Tp = 230 nK. The low- 
est value reached is T/Tp = 0.17(1) with 60 000 atoms 
remaining. Further evaporation does not reduce T/Tp. 

In the case of a single-spin sample of 2.5 x 10^ fermionic 
atoms, we add 6.5 x 10^ atoms of ^^Sr to the dipole trap. 
Evaporation is performed in one single exponential ramp 
over 8 s and results in a pure BEG of 7 x 10^ atoms. To- 
wards the end of evaporation, the fermionic cloud con- 
tains 10^ atoms. The sample is not well thermalized in 
the axial direction of the trap, which we take into ac- 
count by fitting the vertical and horizontal directions 
independently to obtain T/Tp = 0.12 and T/Tp = 0.23, 
respectively. 

Gomparing the ten-state spin mixture to previously 
published results [49, 50 , we reach a ~ 2.5 times lower 
value of T/Tp at a comparable atom number per spin 
state. A deeply degenerate gas is a crucial requirement 
for the reliable operation of a quantum computer based 
on alkaline-earth-metal-like atoms or the exploration of 
SU(A/') magnetism. 



ing thermalization of the sample. In this case, we add 
^^Sr atoms to the system to facilitate sympathetic cool- 
ing. 

In all of these experiments, we determine the atom 
number, temperature T, and Fermi temperature Tp by 
fitting two-dimensional Fermi-Dirac distributions to ab- 



IX. QUANTUM DEGENERATE BOSE-FERMI 
MIXTURES OF SR 

Double-degenerate Bose-Fermi mixtures of ultracold 
gases have been realized with alkali-metal mixtures ^lOTI - 
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Figure 10: (Color online) Double-degenerate Bose- Fermi 
mixtures. The two combinations shown here are ^^Sr + ^^Sr 
(a) and ^^Sr + ^^Sr (b), where all ten nuclear spin states of 
^^Sr are populated. The absorption images of the bosonic 
(left column) and fermionic isotopes (middle column) are 
taken after 25.4 ms of free expansion. The color and length 
scales are identical for all four images. Azimuthally inte- 
grated optical density distributions of the fermionic cloud are 
shown in the right column, together with Fermi-Dirac (solid 
red line) and Gaussian fits (dashed blue line). 



1106] , Li + Yb mixtures jlQ7[ llQ8j , and isotopic mixtures 
of Yb [H], Sr [50], and metastable He [109]. Here, we 
introduce three new quantum degenerate mixtures. We 
combine each of the three bosonic Sr isotopes with a 
ten-spin-state mixture of ^^Sr. The distinguishing prop- 
erty of the three mixtures are the different inter- and 
intraspecies scattering lengths. The interspecies scatter- 
ing length of ^^Sr, ^^Sr, and ^^Sr with ^^Sr are -57 ao, 
162 ao, and 55 ao, respectively, and all allow efficient in- 
terspecies thermalization. 

In a first experiment, we prepare the mixture of 
^^Sr + ^^Sr. In contrast to the experiments presented 
in Ref. [50 and Sec. Villi which used spin-polarized 



fermions, we are now working with the entire ten-state 
mixture of ^^Sr. About 6 x 10^ fermions at 1.15 /jK and 
2.1 X 10^ bosons at 1.1 /iK are loaded into the dipole 
trap. The first of two exponential evaporation ramps 
takes 12 s, during which the onset of BEG in ^^Sr is ob- 
served after 8 s. The BEG is essentially pure after 12 s, 
containing about 3 x 10^ atoms. The fermions remain in 
perfect thermal equilibrium with the bosons. At the end 
of this first ramp, we obtain 7 x 10^ fermions per spin 
state at a temperature of about 100 nK, but still well 
outside the degenerate regime. For further cooling, we 
add a very slow second evaporation ramp of 7.5 s dura- 
tion. Such a slow ramp is required because Pauli block- 
ing decreases the scattering rate between the fermions, 
and superfluidity of the BEG decreases the scattering 
rate between fermions and bosons. The final trap fre- 
quencies are fx = fy = 30 Hz, and fz = 220 Hz. The 
degree of Fermi degeneracy increases substantially to 



T/Tp = 0.15(1) with 15 000 atoms in each spin compo- 
nent at the end of evaporation. The BEG atom number 
decreases slightly to 2 x 10^ atoms [see Fig. 10 'a)]. This 
experiment reaches an eleven-fold degeneracy of distin- 
guishable particles, possibly the largest number of over- 
lapping degenerate gases ever reported. It seems that the 
presence of the bosons reduces the evaporation efficiency 
of the fermions, as deeper degeneracies were achieved in 



absence of the bosons; see Sec. |VIII 

In a second experiment, we use the bosonic isotope 
^^Sr. The bosonic intraspecies and the interspecies scat- 
tering lengths are much larger than in the previous case, 
which we account for by decreasing the density of the 
sample. We load less atoms (1.6 x 10^ atoms of ^^Sr 
at 1.1 /iK and 3.7 x 10^ atoms of ^^Sr at 1.0 /iK) and 
we keep the average trap frequency low by reducing the 
horizontal confinement. We maintain the concept of two 
sequential evaporation ramps of different time constants 
and evaporate slightly deeper than is the previous case, 
but reduce the total evaporation time to 3.4 s. Final trap 
frequencies are fx = 20 Hz, fy = 10 Hz, and fz = 200 Hz. 
This isotopic combination performs worse than the pre- 
vious one, yielding a bosonic sample with only 15% con- 
densate fraction. The horizontal trap frequencies, re- 
quired to be small to keep three-body loss of ^^Sr low, 
does not ensure thermalization of the fermionic sample 
in this direction. We obtain T/Tp = 0.15(5) for the ver- 
tical direction. The BEG contains 5000 atoms, and each 
Fermi sea contains 10 000 atoms. 

As a last experiment, we use the ^^Sr isotope as the 
boson. Starting out with 1.2 x 10^ bosons and 6.5 x 10^ 
fermions both at 1.2 /iK, we reduce the trap depth in 
two ramps of 12 s and 8 s. The atom number of the ^^Sr 
BEG is limited by the negative scattering length, and 
evaporation to a low trap depth is required to remove 
the thermal fraction. We finally obtain a pure BEG of 
4000 atoms immersed in ten Fermi seas, each comprising 
10 000 atoms at T/Tp = 0.11(1) [see Fig.fToFb)]. 



X. CONCLUSIONS 

We have presented a collection of experiments related 
to the creation of BEGs and degenerate Fermi gases of Sr. 
We have shown how electronic and collisional properties 
add to create a strong framework delivering a diversity of 
quantum degenerate samples [see Tab. HI . Key parame- 
ters of Sr quantum gases, such as the BEG atom number 
and the degeneracy of fermionic samples, have been im- 
proved significantly over previous experiments. In par- 
ticular, we have presented the largest BEGs ever created 
by evaporative cooling in an optical trap, we were able to 
reduce the experimental cycle time for BEG preparation 
to 2 s, and we have prepared degenerate Fermi gases at 
T/Tf = 0.10(1). 

Four key properties of Sr make these results possible. 
First, the 461-nm broad-linewidth transition enables effi- 
cient laser cooling. Second, atoms in the metastable 
state can be accumulated in a magnetic trap. This prop- 
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Table III: Overview of the Sr quantum gases and quantum gas mixtures discussed in this article. The isotopes used in a certain 
experiment are designated as "isotope 1" and "isotope 2". If the fermionic ^^Sr isotope is used, the number of populated 
nuclear spin states is given. The initial conditions of the dipole trap are described by the 1/e^-radius of the vertical dipole 
trap beam in the plane of the horizontal dipole trap Wvert and the power of the horizontal and vertical dipole trap beams 
^(horiz/vert), start- The waists of the horizontal dipole trap beam are the same for all experiments, 250(50) /im in the x-direction 
and 18(2) /am in the z-direction. The dipole trap depth is U = ks x i^horiz, start x 5.0(5) /iK/W. The reservoir loading time of 
isotope i is given by tioad,2- After transfer of atoms from the red MOT into the dipole trap, the sample contains Ni atoms at 
temperature T^ of isotope i. Evaporative cooling is performed using one or two nearly exponential evaporative cooling ramps 
with duration tevap,(i/2)- Two ramps are used to create degenerate Fermi gases (DFGs). The second ramp has a much longer 
time constant than the first ramp, which accounts for the reduction of evaporation efficiency by Pauli blocking for samples in 
the quantum degenerate regime. The powers of the dipole trap beams at the end of evaporation are given by P(horiz/vert),end- 
At this point the trap is well approximated by a harmonic potential with trap oscillation frequencies fx,y,z- The last lines of 
the table characterize the resulting quantum gases, their atom number N, the phase transition temperature Tc for bosonic 
isotopes, T/Tf and Fermi temperature Tf for the fermionic isotope, and the peak densities of the gas no. All errors are 
estimated to be below 20 % unless stated otherwise. 
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erty allows us to overcome the low natural abundance of 
^"^Sr. Third, the 689-nm intercombination line, which 
has a linewidth of only 7.4 kHz, enables laser cooling to 
the exceptionally high phase-space density of 0.1. And 
finally, the collision properties of ^^Sr and ^^Sr are favor- 
able for efficient evaporative cooling. Since our first re- 
sults on quantum degenerate Sr, we have implemented a 
range of technical improvements to our apparatus, which 
allow us to make best use of these properties. Besides 
a reduced linewidth of the intercombination line cooling 
laser source and an overall optimization of the exper- 



imental sequence, the most important improvement is 
the design of our dipole trap. Its flexibility allows us to 
quickly choose optimum conditions for each Sr isotope or 
isotope mixture. Some of these techniques might also be 
applicable to other atoms. The metastable state in Yb 
could be used to accumulate atoms, but would require a 
dedicated laser pumping the atoms from the MOT cycle 
into the metastable state. 

The work presented here did pave the ground for the 
creation of bosonic and fermionic Sr lattice gases [58] 
and the creation of Sr2 molecules We plan to use 
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the large ^^Sr BEC as thermal bath for the creation of 
quantum degenerate Rb/Sr mixtures. The exploration 
of S\J{N) magnetism with Sr will greatly benefit from 
the highly degenerate fermionic quantum gases presented 
here. 
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